Stress transmission through three-dimensional granular crystals with stacking faults 
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We explore the effect of stacking fault defects on the transmission of forces in three-dimensional face- 
centered-cubic granular crystals. An external force is applied to a small area at the top surface of a 
crystalline packing of granular beads containing one or two stacking faults at various depths. The 
response forces at the bottom surface are measured and found to correspond to predictions based on 
vector force balance within the geometry of the defects. We identify the elementary stacking fault 
as a boundary between two pure face-centered-cubic crystals with different stacking orders. Other 
stacking faults produce response force patterns that can be viewed as resulting from repetitions of 
this basic defect. As the number of stacking faults increases, the intensity pattern evolves toward 
that of an hexagonal-close-packed crystal. This leads to the conclusion that the force pattern of that 
crystal structure crystal can be viewed as the extreme limit of a face-centered-cubic crystal with a 
stacking fault at every layer. 

PACS numbers: 81.05.Rm, 83.80.Fg, 45.70.Cc 



I. INTRODUCTION 

Static packings of granular particles support stress in 
a nontrivial manner [1]. The forces on any given particle 
from all neighboring particles and from gravity must bal- 
ance vectorally. Stress on a random granular pack is sup- 
ported by a disordered network of these interparticle con- 
tacts in local force balance with each other [2,3]. It has 
been suggested that elastic theory is capable of account- 
ing for these inhomogenaities [4,5]. However, the inde- 
terminacy when the particles are infinitely hard, so that 
particle distortions are not allowed, has led to a search 
for alternative descriptions of the mechanical properties 
of granular materials in the hard-sphere limit [6-9]. 

The proposed alternatives have had varying degrees of 
success. Coppersmith et al. proposed a diffusive q-model 
for the transmission of force through granular packings 
[2,10], and a hyperbolic model was created by Cates et al. 
[11]. Experimentally it has been difficult to rule out any 
of these descriptions conclusively, but recently progress 
has been made in studying the conditions under which 
these models apply. 

Much understanding and additional surprises have 
been gained from experiments studying the response to 
a localized external force, similar to a Green's function. 
DaSilva and Rajchenbach [12] studied a two dimensional 
packing of rectangular photoelastic bricks. When ap- 
plying an external force at the top of the packing they 
measured the response forces between neighboring bricks 
within the pack. The results showed that at a given 
depth the maximum interparticle force was located di- 
rectly below the position of the applied force. Deeper in 
the packing (farther from the applied force) they found 
the width of this peak of maximum force to broaden as 
the square root of the depth. Similar results were found 
by Moukarzel et al. in recent displacement response ex- 



periments [13]. Such findings support diffusive models. 
However, a cross-over to a linear widening of the response 
peak cannot be ruled out for deeper packings. Reydellet 
and Clement looked at the response forces at the bot- 
tom surface of an amorphous three-dimensional packing 
of spheres [14]. Again they found the maximum force 
situated directly beneath the applied force, but now the 
width of the peak grew linearly with depth from the sur- 
face. This linear growth lends credence to the elastic 
theories. Similar elastic-type behavior was found in two- 
dimensional packings by Ceng et al.. In addition those 
researchers also found a strong dependence on the spatial 
ordering of the particles [15,16]. 

The importance of spatial ordering was further shown 
in three dimensional granular crystals in Ref. [17]. 
Within face-centered-cubic (FCC) crystals, the force was 
found to be supported along straight lines of contacts be- 
tween beads, resulting in three areas of large force at the 
bottom surface in response to a locally applied external 
load at the top surface. In contrast, an hexagonal-close- 
packed (HCP) arrangement resulted in a ring of large 
force at the bottom surface. In both cases there was a 
local minimum of force directly beneath the point of ap- 
plication. Such patterns are reminiscent of the hyperbolic 
models, which predict straight-line propagation of forces. 
Bouchaud et al. have extended these hyperbolic models 
to include the splitting of force chains at defects within 
a granular packing [18]. Such splittings are predicted to 
lead to the elastic-like behavior that has been observed 
in amorphous packings. 

It is, therefore, of interest to explore the regime be- 
tween perfectly-ordered crystalline bead packs and com- 
pletely amorphous arrangements. As a first step in this 
direction, we present an experimental study of the ef- 
fect on the force propagation of a small number of stack- 
ing fault defects within otherwise perfectly ordered FCC 
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crystals. 

II. EXPERIMENTAL METHODS 

Employing the triangular acrylic cell described by Blair 
et al. [20] and Mueggenburg et al. [17] and a 'by hand' 
construction technique, we produced large FCC crystals 
of approximately 20, 000 soda lime glass spheres of di- 
ameter 3.06 ± 0.04mm with an extremely low number 
of defects (estimated to be fewer than ten beads signifi- 
cantly out of place in the entire pack) oriented as horizon- 
tal planes of triangular order. Furthermore, we carefully 
controlled the stacking of the planes in order to create 
true three-dimensional crystals and to introduce stack- 
ing faults at specific locations. 

In this manner, we constructed five different granular 
FCC crystals, each having one or two defects in the stack- 
ing order. We then applied an impulsive force to a small 
region, approximately two beads in diameter, centered on 
the top of the crystal [17]. Using a carbon paper tech- 
nique, we characterized the normal forces on the top and 
bottom surfaces of the crystal in response to this force 
[2,3,17,19,20]. A piece of carbon paper was sandwiched 
between the crystal and a piece of white paper on both 
the top and bottom surfaces of the crystal. The applied 
force caused individual beads to press into the carbon 
paper and leave marks on the white paper. These marks 
were then digitized and image analysis software was used 
to calculate the area and intensity of each mark [17]. 

We conducted twenty experimental runs on each crys- 
tal. For each run the resultant carbon marks were aligned 
with the expected triangular lattice. The areas and inten- 
sities were averaged over experimental runs and, in order 
to improve statistics, were averaged over symmetries of 
the crystal (one refiection and two rotations) when such 
averaging did not qualitatively alter the intensitiy pat- 
terns. 



III. RESULTS 

A pure FCC crystal consists of layers of triangu- 
lar order stacked in such a way that every third layer 
lies on top of the first. We annotate this stacking 
as ABC ABC ABC ABC ABC. . . A stacking fault corre- 
sponds to a disruption of this order and may be consid- 
ered to be a boundary between two pure FCC crystals 
with different stacking orders. 

We constructed a 19 layer crystal with a stacking fault 
centered at the tenth layer. This stacking is represented 
as ABCABCABCA CBACBACBA The underline rep- 
resents one pure FCC crystal, while the overline corre- 
sponds to a second pure FCC crystal with a different ori- 
entation. At the boundary between the two crystals there 
exists an A layer which correctly matches both stacking 




FIG. 1. 19 layer FCC crystal with one stacking fault at 
layer 10. a) Schematic representation of layer stacking show- 
ing a boundary between two FCC crystals, b) Force intensity 
plot at bottom surface in response to a quick impulse applied 
to a small region at the center of the top of the packing, c) 
Force intensity plot after averaging over crystal symmetries, 
which accentuates anomalous interior features. 

orders and can thus be considered to be in both crys- 
tals. Note that, by necessity, the three layers centered 
around the stacking fault CA C follow an hexgonal-close- 
packed stacking order. When applying a quick impulse 
to a small area at the center of the top of this crystal, the 
bottom surface displayed six regions of large force at the 
vertices of a regular hexagon as shown in figure 1. The 
regions of large force are approximately the same size 
as the area over which the impulse was applied on the 
top surface. We observed another region of larger than 
background force inside the hexagon, but closer to one 
face. At other symmetric places of the pattern, areas of 
larger than background force are present, but with much 
lower intensity. It is not clear if this area of heightened 
force is a real feature of the crystal with stacking fault or 
if this is merely an extreme fiuctuation. Averaging over 
crystal symmetries would accentuate such fiuctuations, 
and thus we also show the intensity pattern without this 
averaging. 

Additionally, we explored possible assemblies of crys- 
tals with two stacking faults. Figure 2 represents a 13 
layer FCC crystal with stacking faults at the seventh 
and eleventh layers. The stacking order (from top to 
bottom) was ABC ABC A C B A C AB as shown schemati- 
cally in figure 2(a). We found that this crystal produced a 
twelve-peaked intensity plot (Figure 2b). The plot shows 
two pairs of peaks of nearly equal height aligned along 
the faces of an equilateral triangle. Each pair of peaks 
is approximately four beads wide and the separation be- 
tween pairs is approximately three beads. 

Another crystal with two stacking faults consisted of 
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FIG. 2. a) Schematic representation of layer stacking of 
a 13 layer crystal with faults centered at layers seven and 
eleven, b) Intensity plot of forces on the bottom surface show- 
ing twelve peaks of high force. 



FIG. 4. a) Representation of layer stacking for a 13 layer 
crystal with stacking faults at the fifth and ninth layers, b) 
Intensity plot of forces on the bottom surface showing three 
larger peaks each flanked by a pair of smaller peaks. 




FIG. 3. a) Schematic representation of layer stacking of an 
11 layer crystal with faults centered at layers four and eight, 
b) Intensity plot of forces on the bottom surface showing three 
large peaks each flanked by a pair of smaller peaks. 



eleven layers with faults centered at the fourth and eighth 
layers ( ABC A C BA G ABC ) . For this crystal, we found 
three tall peaks flanked by peaks half the height of the 
taller ones arranged along the faces of an equilateral tri- 
angle (Figure 3). The width of each of the taller peaks 
was slightly larger than that of the smaller flanking peaks 
which correspond to the width of the area over which the 
force was applied at the top surface. 

The next crystal we built contained thirteen layers 
with faults at the fifth and ninth layers. The stacking 
order was ABC AB AC B ABC AB . We found an inten- 
sity pattern (Figure 4) similar to that of the previous 
crystal with two faults. This crystal produced three tall 
peaks flanked by peaks about one third the height of the 
taller peaks arranged along the faces of an equilateral tri- 
angle. The taller peaks were approximately three beads 
wide and the shorter peaks were approximately two beads 
wide. In contrast to figure 3, which corresponds to a crys- 
tal with fewer layers above and below the stacking faults, 
the flanking peaks are now farther away from the taller 
peaks. 

Another possible type of stacking fault involves the 
omission of one plane from a large FCC crystal. This 
fault maintains the orientation of the crystal above and 
below the fault and thus might be considered a more basic 
defect than the grain boundary described above. How- 
ever, this fault could be described as a special case of two 
grain-boundary faults. Figure 5 shows the results for a 
crystal with layer stacking ABCABC'ACABCABC . Al- 




FIG. 5. a) Layer stacking for a 14 layer crystal with one 
missing layer fault - also described as two grain boundary 
stacking faults at layers seven and eight, b) Intensity plot of 
forces on the bottom surface showing six peaks in a triangular 
arrangement. 



though one might suppose that this crystal exhibited the 
least amount of stacking disorder, we found that the force 
intensity pattern consisted of six approximately equal- 
height peaks, one at each of the vertices of an equilateral 
triangle and one at each of the midpoints of the sides 
(Figure 5). The peaks at the vertices of the triangular 
pattern are four beads in width while those on the sides 
are three beads wide. We describe this pattern below as 
a special case of the nine or twelve peaks shown in the 
crystal packings with two faults. 



IV. DISCUSSION 

As in the experiments performed on pure FCC crys- 
tals by Mueggenburg et at. [17], we can reconstruct the 
force chain paths through the crystal based on the force 
intensity plots from the experiments. More importantly, 
the simple force balance model used to explain the force 
patterns in that study, also explains the present results. 

Forces will travel through pure FCC sections of the 
crystal, along the spines of a tetrahedron. Thus, to un- 
derstand an elementary stacking fault, all we need to 
determine is what happens at the interface between two 
FCC crystals with different orientations. Again, the force 
balance model can be applied. In order to maintain force 
balance at the interface, the force chain must split. Note 
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that locally this is the same geometry as an HCP crys- 
tal. Below the stacking fault each of the new forces will 
travel through the lower portion of the crystal in straight 
line paths. Figure 6(a) illustrates this force splitting at 
a grain boundary. 

We found a regular hexagonal pattern of force in our 
first crystal because the fault was placed in the middle 
of the packing. In general, we would expect a six-peaked 
pattern in which two peaks reside along each side of a 
triangle. As the defect is moved higher or lower in the 
crystal, these pairs of peaks would move towards or away 
from each other. 

It is simple to extend this idea to crystals contain- 
ing two stacking faults. The crystal can be thought of 
as three pure FCC crystals with top and bottom sec- 
tions having the same stacking order and the middle sec- 
tion having a different stacking order. At each interface 
the force splitting process is repeated resulting in twelve 
peaks in the intensity pattern (3 split into 6 and then 
split into 12). Figure 6(b) illustrates this multiple split- 
ting process. However, most of our results do not demon- 
strate the expected twelve-peaked pattern. This can be 
explained by force recombination. 

For certain combinations of stacking faults, the forces 
will converge upon the same lattice point on the bot- 
tom of the crystal. In other cases, two forces could 
travel along parallel lines through the crystal and ar- 
rive at neighboring lattice points on the bottom surface. 
In the former scenario, the force intensity at the con- 
vergence point should be twice the intensity of a lattice 
point with only one line of force. In the latter scenario, 
the forces should produce a wider than normal peak (the 
exact width would depend upon the width of the applied 
impulse and the separation of the force chains within 
the packing). These two effects account for our observed 
force patterns from crystals with two faults. 

In the second and third crystals with two stacking 
faults (Figures 3 and 4), the taller central peaks are places 
where forces have recombined and the smaller flanking 
peaks correspond to single force chains. The widths of 
the peaks corroborate this explanation in both cases. The 
fourth example, where one plane of the lattice has been 
removed, can be categorized as a crystal with two stack- 
ing faults in which the twelve peaks have been reduced 
to six due to this recombination. The corner peaks in 
the intensity plot (Figure 5) are caused by parallel forces 
terminating at adjacent lattice points and the side peaks 
are caused by two forces terminating at the same lattice 
point as shown in figure 6(c). The widths of the peaks 
support this explanation. This is the extreme case of the 
trend from figure 3 to figure 4 where the flanking peaks 
moved towards the corners of the triangular pattern as 
the spacing between defects relative to the overall crystal 
height decreased. In the case, where the two defects oc- 
cur at successive layers, the flanking peaks are positioned 
very near the corners and thus overlap with each other. 




neighboring points 

FIG. 6. Diagrams of force transmission through a crystal 
packing with stacking faults, a) Forces must split at the stack- 
ing fault, b) Multiple stacking faults cause multiple force 
splittings, c) In some cases multiple force chains may recom- 
bine or otherwise arrive at the bottom surface at neighboring 
points and thus produce fewer areas of large force. 
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V. CONCLUSIONS 

In this study, we have estabhshed how stacking faults 
affect the transmission of stress through an FCC crys- 
tahine granular assembly. We have observed that the 
force chains split into two equal components at the 
boundary between different FCC crystal orientations. 
However, we note that the forces are transmitted as in 
a pure FCC crystal before and after splitting. Our find- 
ings reinforce those in reference [17] and the force balance 
model proposed therein. 

We have identified the elementary stacking fault de- 
fect as the boundary between two pure FCC crystals of 
different stacking orders. All further stacking faults can 
be described as repetitions of this basic defect. Alter- 
nately, we can express the stacking faults as sections of 
HCP crystals within the overall FCC crystal. Although 
we find the grain-boundary description the most lucid, 
this other interpretation is geometrically equivalent. 

As the number of stacking faults increases, the inten- 
sity pattern looks increasingly like that of an HCP crys- 
tal, and indeed the local stacking order around a fault is 
HCP. This leads to the conclusion that the force pattern 
of an HCP crystal can be viewed as the extreme limit of 
an FCC crystal with a stacking fault at every layer. 
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